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There is a sharp increase in the relative number of degrees 
of freedom at the deconfining phase transition. Characteriz- 
ing this increase using the Polyakov Loop model, we find that 
for a nearly second order deconfining phase transition, the 
medium-induced energy loss turns on rapidly above T c , pro- 
portional to the relative number of degrees of freedom. Fur- 
ther, energy loss is logarithmically dependent on the screening 
mass, and thus is sensitive to nearly critical scattering. 

Experiments indicate that for the central collisions of 
large nuclei, A ~ 200, there are marked changes between 
energies of yfi/A = 17 GeV, at the SPS, and 130 GeV, 
at RHIC ffl . Comparing central AA collisions to pp, the 
spectrum of semi-hard particles is rather different. At the 
SPS, in AA the hard pt spectrum, scaled by the number 
of binary collisions, is enhanced over pp. At RHIC, the 
opposite is true: the semi-hard p t spectrum per nucleon- 
nucleon collision, is suppressed in central AA, relative 
either to peripheral AA, or pp 0]. This could be the 
result of "energy loss" where a fast colored field 

loses energy as it passes through a thermal bath. In pe- 
ripheral A A collisions, secondary hadrons are distributed 
anisotropically in the transverse momentum p t M . Ex- 
perimentally, this azimuthal anisotropy increases with p t 
until pt ~ 2 GeV, at which point it flattens 0. This 
flattening may also be due to energy loss || . 

In the limit of infinitely large nuclei, A — > oo, it is 
plausible that the initial energy density produced in a 
central A A collision — at a fixed value of y/s/A — evolves 
into a system in equilibrium at a temperature T. With 
great optimism, assuming that A ~ 200 is near A = oo, 
one might imagine that the difference between SPS and 
RHIC is because temperatures reached at RHIC exceed 
T c , the critical temperature for QCD. 

Thus it is of interest to know how quantities change 
as one goes through the phase transition. In this paper 
we give an analysis in terms of the Polyakov Loop model 

In QCD, there is a large increase in the number of 
degrees of freedom at the deconfining phase transition. 
We count degrees of freedom as appropriate for the pres- 
sure of free, massless fields at nonzero temperature, so 
if each boson counts as one, then each fermion counts 
as 7/8. In the hadronic phase, pions contribute = 3 
ideal degrees of freedom. By asymptotic freedom, at infi- 
nite temperature QCD with three flavors of quarks is an 
Ideal Quark-Gluon Plasma, with cqgp = 47 \ degrees of 



freedom. This is an increase of more than a factor of ten. 

To measure the change in the number of degrees of 
freedom, we introduce the relative pressure, n(T): at 
a temperature T, this is the ratio of the true pressure, 
p(T), to that of an Ideal Quark-Gluon Plasma, Pidcai = 
cqgp(7tV90)T 4 : 
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By asymptotic freedom, QCD is an ideal gas at infinite 
temperature, and so 



n(oo) = 1 



(2) 



For T < oo, corrections to ideality are determined by the 
QCD coupling constant, a s oc 1/ log(T), with n(T) — 1 oc 
-a. @. 

For an exact chiral symmetry which is spontaneously 
broken by the vacuum, about zero temperature the free 
energy is that of free, massless pions. Thus at zero tem- 
perature, the relative pressure is the ratio of the ideal gas 
coefficients 1131]: 
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cqgp 
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At low temperature, corrections to ideality are given by 
chiral perturbation theory for massless pions, n(T) — 
n(0) ~ +(T/f 7T ) 4 n(0), with / w the pion decay con- 
stant. In QCD, pions are massive, and the relative 
pressure is Boltzmann suppressed at low temperature, 
n(T) ~ exp(-m 7r /T)(m w /T) 5 / 2 , so n(Q) = 0. 

Given the great disparity between c n and cqgp, con- 
sider an approximation where the hadronic degrees of 
freedom are neglected relative to those of the deconfmed 
phase [ fl4| ]. Then the relative pressure vanishes through- 
out the hadronic phase, n(T) = for T < T c . The 
question is then: how does the relative pressure go from 
zero at T c , when deconfinement occurs, to near one at 
higher T? 

This can be answered by numerical simulations of Lat- 
tice QCD []l5]]. Consider first quenched QCD, with pure 
glue and no dynamical quarks, which is close to the con- 
tinuum limit fusfl . For three colors, the Lattice finds no 
measurable pressure in the hadronic phase (glueballs are 
heavy), so our approximation of n(T) = when T < T c 
is good. n(T) increases quickly above T c , and is ~ .8 
by T ~ 2T C . To characterize the change in the relative 



1 



pressure, consider the ratio of e — 3p, where e(T) is the 
true energy density of QCD, to the energy of an Ideal 
Quark-Gluon Plasma, eideai = 3p;deai : 



m e (T)/T ot£ Q ~ n 1/4 



(8) 



3p 



e idcal 



T dn 
3df 



(4) 



Lattice simulations find that this ratio has a sharp 
"bump" at ~ 1.1T C , suggesting that the relative pressure 
changes quickly, when the reduced temperature, 



t 



(5) 



is small, t ~ .1. 

The Lattice is more uncertain with dynamical quarks. 
The pions are too heavy, and it is not near the continuum 
limit. So far, the Lattice finds that n(T/T c ) is about the 
same with dynamical quarks as without p5 17 1. This 



suggests that the pure glue theory may be a reasonable 
guide to how the relative pressure increases above T c . 
The approximate universality of n(T/T c ) is remarkable. 
At present, the Lattice finds no true phase transition 
in QCD, with T c smaller by ~ .6 than in the quenched 
theory Indeed, even the ideal gas coefficients are 

very different: cqgp is only 16 in the quenched theory, 
versus 47 \ in QCD. 

The greatest change with dynamical quarks is a small, 
but measurable, pressure in the hadronic phase. While in 
the quenched theory n(T) ~ for T < T c , with dynami- 
cal quarks, although n(0) ~ 0, there is a nonzero relative 
pressure at the critical temperature, with n(T c ) ~ .1 [15[ . 
Indeed, with no true phase transition, an approximate T c 
can only be defined as the point where the relative pres- 
sure increases sharply, reaching n ~ .8 by 2T C |jl5| . 

The Polyakov Loop model [p|-pT[ is a mean field the- 
ory for the relative pressure. In a pure glue theory, the 
expectation value of the Polyakov Loop, lo(T), behaves 
like the relative pressure: it vanishes when T < T c , and is 
nonzero above T c . Indeed, again by asymptotic freedom, 
£o — > 1 as T — > oo. The simplest guess for a potential for 
the Polyakov Loop is: 



v(t) = - b p\ 2 + ~ 



(6) 



Defining £o as the minimum of V(£) for a given hi (T), 
the relative pressure is given by p-O]: 



n(T) = -W(£ ) = £ 4 



(7) 



b 2 > above T c (b 2 (T) -> 1 for t -> 00), and < be- 
low T c . Thus if the relative pressure changes when the 
reduced temperature t ~ .1, the change for £q(T) ~ n 1 / 4 
is even more rapid, within 2.5% of T c . 

For two colors, (||) is a mean field theory for a sec- 
ond order deconfining transition ]T^ |. The £ field is 
real, and so the potential defines a mass: (mi/T) 2 = 
(l/Z s )d 2 V/d£ 2 , with 



where Z s is the wave function normalization constant for 
£, Z s = 3/g 2 , up to corrections of order g° This is 

measured from the two point function of Polyakov loops 
in coordinate space, oc (1/r) exp(— mgr) as r — > 00. 

For three colors, £ is a complex valued field, and a term 
cubic in £ appears in V(£), —b3(£ 3 +£**)/ 6. This produces 
a first order deconfining transition, where £q jumps from 
at T~ to £ c = 263/3 at T+ |o|. The £ field has two 
masses, from its real (mi) and imaginary (frig) parts. At 
^/~Z~ s mi/T = £ c ; from the Lattice, ^J~Z~ s mt/T ~ .3 
|15| , which gives 63 ~ .45. This small value of 63 re- 
flects the weakly first order deconfining transition for 
three colors EJl(| . The mass for the imaginary part 
of £ is y/Z~mz(T)/T oc y/b^£ - n 1 / 8 ; at T+, m e / mi = 3. 
With dynamical quarks, in principle a term linear in £, 
-b x (l + e*)/2, can also appear in V(£) (2f|. If the pion 
pressure is included below T c , however, b\ is very small, 
< .03. 

Thinking of £ provides a useful way of viewing the 
deconfining phase transition. For a strongly first order 
transition — as appears to occur for four or more colors 
fHf — £q, jumps from zero below T c , to a value near one 
just above T c . As £0 is near one, the deconfined phase is 
presumably well described as a nearly Ideal Quark-Gluon 
Plasma |2^] . In this case, there is a hadronic phase below 
T c , and a Quark-Gluon Plasma from T c immediately on 
up. 

In contrast, for three colors the deconfining transition 
is weakly first order. As the energy density is discontin- 
uous at T c , for small t the relative pressure is linear in 
the reduced temperature, 



n(T) ~ 3rt ; 



(9) 



here r = e(T+)/eideai(7c) is the ratio of the energies at 
T c , in the deconfined phase versus an Ideal Quark-Gluon 
Plasma. For quenched QCD, r ~ 1/3 |l6|], which gives 
n(T) ~ t, and so £q(T) oc i 1 / 4 . Except very near T c , this 
simple estimate agrees with more complicated analysis 
using 63 =/= ) . For example, at only 5% above T c , 
this estimate gives £q ~ .05 1 / 4 ~ .5. For three colors, 
then, there is a (non)-Ideal Quark-Gluon Plasma only 
at temperatures above ~ 2T C ; between T c and ~ 2T C , 
the Polyakov Loop dominates the free energy, going from 
- .5 at 1.05T C to - 1 by 2T C . 

The difference between these two scenarios: a strongly 
first order transition, where £o(T) is approximately con- 
stant above T c , and nearly second order behavior, where 
£q(T) changes significantly, is in principle observable. As 
an example, consider energy loss for a fast parton, with 
a high energy E. We first give a general discussion of 
energy loss in a medium and then discuss the dif- 

ferences between a strong first order transition, and one 
which is nearly second order. 

We introduce the energy scale || , 

E CI = ^L 2 , (10) 
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where A is the mean free path and L is the thickness 
of the medium. The high-energy jet loses energy by ra- 
diating gluons with energy to < E. There are several 
contributions to the total energy loss of the jet, AE, de- 
pending on the energy of the radiation. For the contri- 
bution from lo > E'er, which exists if E > E cr , effectively 
only one single scattering occurs (this is the so-called 
factorization regime) and so that contribution is medium 
independent In what follows we rather focus on the 
medium-induced energy loss, from the region where u> is 
less than E CI . 

For very small frequency, lu < -Elpm = Xmj, the 
formation time t/ ~ Lu/mf of the radiation from 

the hard jet is short, and so incoherent radiation takes 
place. This is the so-called Bethe-Heitler regime; the 
contribution to AE is just a sum from single scatter- 
ings on L/X scattering centers. In the high-energy limit 
E, E CI ^> Slpm the region of phase space with u> < E^pm 
contributes little to AE and will be neglected. 

The largest contribution is rather from the Landau- 
Pomeranchuk-Migdal (LPM) regime, where successive 
scatterings coherently interfere ■ Integrating the ra- 
diation intensity distribution over lu from zero to some 
energy E* yields a total energy loss of G|j5| 



AE 



3 a. 



VE~E* log 



2E 



(11) 



There is a logarithmic sensitivity to the infrared scale 
mi Qj2^] . When the jet energy E is less than the factor- 
ization scale E cr , we can integrate lu all the way up to 
E* = E, so 



3a s > 2E 

AE a/ E E cr log — — 

71 lj Jll i) 



(E < E cr ) . (12) 



Note that — AE should not exceed E. This requires that 
me is not so small that the logarithm overwhelms the 
~ a s . However, E cr is small near T c , so minijets with 
energies of at least a few GeV are above E cr anyways, 
and eq. ( |l2| ) does not apply. 

Rather, for jet energies greater than E CT , the total 
medium-induced energy loss is given by integrating over 
u> up to the factorization scale E CT ] setting E* = E CI 
in flub, 



3a s 2E 

-AE '- E CT log — - 

7r Lmi 



(E > E CI ) . (13) 



To compute the critical energy E CI , we need the inverse 
mean free path, A^ 1 . This is approximately the product 
of the density, p, times the elastic cross section, a c \. The 
elastic cross section is quadratically divergent in the in- 
frared. This divergence is naturally cut off by mi, so the 
clastic cross section <r c \ cx a 2 s /mf, and 



oc p 



(14) 



The scale E cr is then proportional to p; this follows au- 
tomatically from our assumption that A -1 ~ po e \- Thus 
in the high-energy regime above E cr , ignoring the loga- 
rithmic dependence upon me, energy loss is proportional 
to p; below that scale, to ^fp. This has been emphasized 
by Baier, Dokshitzer, Mueller, and Schiff ||. 

Before giving estimates of me and A, we can under- 
stand how energy loss changes, depending upon the order 
of the deconfining phase transition. For a nearly second 
order transition, me/T is small near T c , eq. and then 
increases rapidly. As the energy loss AE depends loga- 
rithmically on me, for small me energy loss is enhanced. 
This is directly analogous to critical opalescence. In con- 
trast, for a strongly first order transition, me is large at 
T c + , with me/T approximately constant with increasing 
temperature. 

What are reasonable values for me and A ? In the 
extreme perturbative regime, T S> T c , (static) electric 
fields are heavy, with a mass oc ifa^T, while the static 
magnetic fields are light, m mag oc a s T. The inverse mean 
free path, A -1 , equals the damping rate for a gluon with 
momentum ~ T, and is oc a s T. 

At temperatures ~ 2T C , this ordering is reversed, as 
static electric fields are significantly lighter than static 
magnetic fields: me ~ 2.5T, while the static mass for 
magnetic glueballs is m 



n) ag 



6T |24|. There are no 
estimates of the damping rate; we guess that 7 ~ A -1 ~ 
T. This seems reasonable for a quasiparticle with such 
a mass, as j/me ~ 1/2.5 = .4 is less than one. If the 
width were much larger, then it would not make sense to 
speak of quasiparticles. Conversely, in a strongly coupled 
system, it is unreasonable to think that the width could 
be much smaller than the mass. 

In the derivation of energy loss, implicitly it is assumed 
that multiple scatterings of the hard jet are independent 
of each other. This requires that the range of the po- 
tential is smaller than the mean free path, mj l < A 
|^,||, which is equivalent to having quasiparticles with 
relatively narrow width. 

At 2T C , then, mj/\ - (2.5T) 2 T. Below 2T C , by our 
mean field analysis, then, mf/X ~ 6.25T 3 n(T). Notice 
that as mf ~ \fnT 2 , and 7 = A -1 ~ \/nT, that j/me ~ 
n 1 / 4 : the £ quasiparticles become narrower as T — > T c + . 
Consequently, for temperatures near T c , 



E CI - 20.GeV x n{T) 



(15) 



This number is a very crude estimate, obtained by taking 
L ~ 5 fm and T c - 175 MeV. For QCD, n(T c ) - .1, so 
then E cr - 2 GeV. 

It is interesting to note that while the system is be- 
coming increasingly dilute near T c for a nearly second 
order transition, typical minijets with E on the order 
of several GeV automatically have energies larger than 
E CI . This implies that i) some part of their energy loss 
is medium independent, from the factorization regime 
lu > E cr ; and that ii) their medium-induced energy loss is 
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given by eq. (|lq), not (12). Physically, this is because the 
formation time tf ~ ui/mj of the radiation from the hard 
jet grows like n -1 / 2 as T — > T+. In all, -E cr is small near 
T c , and the regime where A_E scales with ,/p shrinks. 
That regime only emerges as E cr increases: by 2T C , with 
n(2T c ) ~ 1, _E cr has risen dramatically, to ~ 160 GeV ! 
Thus, the energy loss ( filf ) turns on very rapidly as T in- 
creases from T c . In part this is because of the factor of 
T 3 in eq. (fL5|), which is present whatever the order of the 
phase transition. For gauge theories with three colors, 
such as QCD, however, additionally there is an increase 
cx n{T) gi{l7). 

With these numbers, the logarithmic sensitivity to the 
changes in the screening mass can be significant. For 
E = 25 GeV, L = 5 fin, T = T c , then if m e ~ 2.5T C , 
\og(2 E / (Lmf)) ~ 2.3. Including the change in the 
screening mass near T c , with mi ~ 2.hT c n 1 / i , the log- 
arithm changes to \og(2E / (Lmj)) ~ 3.5. This is an in- 
crease by about 50%. At smaller E, the sensitivity to me 
is even stronger. 

Even with the limitations of our approximations, it is 
clear that since the density vanishes as T — ► T c + , any con- 
tribution from the deconfined phase vanishes like some 
power of n(T) |p5[ . For example, estimates indicate that 
dilepton production from the (nearly ideal) deconfined 
phase [E6| is about as large as that from the hadronic 



phase 



Near T c , dilepton production from the decon- 



fined phase should be strongly suppressed, quadratically 
in the density: 



dN P 



d 4 a 



2 rr.4 

cx n 1 



(16) 



This assumes that the density scales as the relative pres- 
sure, n(T). Whatever the exact form, it is clear that 
dileptons from the deconfined phase will turn on later 
(in T, or */sJA) than energy loss, which close to T c is 
linear in the density, and thus in n(T). 

The amplitude for dilepton production involves a vir- 
tual photon off the light cone; this has no power like 
infrared divergences at either leading or next to leading 
order in a s , and so dilepton production should not be 
especially sensitive to the screening mass, me. This may 
not be true for the production of real photons; this in- 
volves processes on the light cone, with severe infrared 
divergences at next to leading order |28|. At leading or- 
der in a s , real photons from the deconfined phase will be 
suppressed like ~ n 2 , as in eq. (fuf). It is conceivable that 
infrared effects from higher order diagrams change this to 



n /mf 



,3/2 



If true, then there is a sequential series 



of effects: increasing T from T c , as the relative pressure 
increases, first energy loss turns on, then perhaps real 
photons, then dileptons. 

In closing, we remark that changes in the screening 
mass could also be significant for the collision rates of 
quarks and gluons. This is particularly true for quantities 
like the elastic cross section, the color conductivity, etc, 
which are sensitive to me p9j. 
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